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[^ ■ ABSTRACT 

^ . An analytic collapse model for the formation and density distribution of virialized 

cold dark matter halos is proposed. Hierarchical structure formation is taken into 
account explicitly. Monte Carlo methods are used to generate samples of mass histories 

Qv ■ of virialized halos. The mean density distribution found from the collapse model is in 

0> ■ good agreement with numerical results in the mass range from 10^^ Mq to IO^^Mq and 

^ I in the radial range form 0.05r2oo to r2oo- 
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1. Introduction 

The current standard paradigm of cosmolog- 
ical structure formation is summarized by three 
main points: The Universe on large scales is dom- 
inated by collisionless dark matter; the seeds of 
collapsed systems are small-amplitude Gaussian 
fluctuations of the uniform primordial density 
field and the evolution towards virialized systems 
occurs by gravitational instability. In hierarchi- 
cal models of structure formation, the amplitude 
of the primordial density fluctuations increases 
with decreasing scale. This leads to a picture 
where small objects form first and, by subsequent 
merging, form larger objects later on. Numerical 
simulations (e.g. Lacey & Cole 1994) confirmed 
this picture to a large extent. 

A striking outcome of numerical simulations 
is the 'universal' shape of the profile of virialized 
density distributions e.g. of the form (Navarro et 
al. 1996, NFW96 hereafter) 
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where ro is a characteristic length scale and po a 
characteristic density. Recent numerical simula- 
tions (e.g. Moore et al. 1998) indicate a some- 



what steeper innermost slope p 



r^ f 



-1.4 



How- 



ever, in the radial range [0.05r2oo 1^200] a broad 
agreement is found. 

Early analytical studies of the density distri- 
bution of collapsed objects are based on spher- 
ically symmetric systems (Gunn & Gott 1972, 
Gunn 1977) forming around sharp density peaks. 
Gaussian initial conditions have been taken into 
account by Hoffman & Shaham (1985). These 
studies suggest a virialized density distribution 
of the form p oc r~'^ with 7 ~ 2. From the 
point of view of the current standard paradigm 
of structure formation these models suffer from 
one main defect: They assume continuous and 
self-similar infall. There have been several at- 
tempts to develop more consistent models. Most 
of these studies seek to explain the form of the 
density distribution (||). For example, it has been 



argued (Syer & White 1998) that it arises from 
repeated mergers. Semi-analytical (Avila-Reese 
et al. 1998) and combined analytical and numer- 
ical approaches have been conducted as well. By 
including dynamical friction and tidal stripping, 
Nusser & Sheth (1998) e.g. find that the final 
density profile mainly depends on the rate of in- 
fall. The scaling of the density distribution (Q) 
with formation redshift has been investigated by 
Raig et al. (1998). 

One of the key issues regarding hierarchical 
structure formation is how well the formation his- 
tory is remembered by a collapsed system. The 
aim of this letter is to investigate this question. 
As a working hypothesis we assume that each 
step of the formation history is remembered by 
the conservation of the corresponding total en- 
ergy of infalling matter. Under this assumption, 
a simple spherical symmetric collapse model is 
proposed. The model is developed along the line 
of earlier analytical studies (e.g. Gunn 1977, 
Hoffman &: Shaham 1985). However, the con- 
dition of self-similarity is abandoned and the hi- 
erarchical formation history is take into account 
explicitly. 

2. The Model 

The model is based on a) formation history 
and b) a collapse model describing how infalling 
matter virializes. The cosmological context is im- 
plicitly present in a). 

2.1. Formation History 

The excursion set formalism (Bond et al. 1991, 
Lacey &: Cole 1993, 1994) provides an analytical 
model of the merger history of a collapsed ob- 
ject. Given the mass M^ of an object at red- 
shift Zfc, it predicts the probability distribution 
of the mass M^j^i of the progenitor at an earlier 
epoch Zfc+i- Monte Carlo methods can be used 
to generate particular merger histories. In doing 
so, we follow here the method of Eisenstein & 
Loeb (1996). Starting from a halo mass Mq at 



zq = 0, the merger history is traced by stepping 
backwards in time until a typical mass marking 
the seeding of the object is reached. A particular 
mass history then is described by a set {{rrik, Zk), 
k = 0,1,2, ... , K} where Tiik = M^ — Mk+i is the 
mass falling in at redshift z^ and K denotes the 
total number of formation steps. 

2.2. Total Energy of Infalling Matter 

In the spherical symmetric case, the total en- 
ergy Ek, the infalling mass ruk and its turnaround 
radius Rm are related by 
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where Mk+i is the mass interior of Rm, i.e. the 
mass collapsed at z > Zk- For ruk <C MkJ^-i, Rm 
is approximately 
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where pm is the mean density of the region of 
mass Mfe_|_i at maximum expansion. For an 
Einstein-de Sitter cosmology (0 = 1, A = 0) 
considered here (e.g. Peebles 1980) 
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Here pQ is the background density at the current 
epoch 2 = and Zk+i is the redshift at collapse. 
Prom eqs. (||) and Q the total energy Ek of 
matter of mass nik falling in at formation step k 
is 



Eu 
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How does this approximation of Ek fit the general 
picture of hierarchical structure formation? The 
basic assumption which has been made is that 
the density distribution of collapsed matter is ap- 
proximately spherical symmetric, an approxima- 
tion which is reasonable for virialized systems. 
No assumption about the distribution of the in- 
falling matter has been made, however. The in- 
falling matter can be distributed continuously or 



it can be concentrated in one or several spot at 
r = Rm- In this sense the approximation for Ek 
accounts for a hierarchical structure formation 
process where larger objects form by the collapse 
of already collapsed smaller objects. 

2.3. Collapse Model 

The working hypothesis for the proposed col- 
lapse model is that a virialized halo has a mem- 
ory of its formation history in the sense that 
at each formation step the total energy of in- 
falling matter is conserved. To determine the 
virialized structure, two additional simplifying 
assumptions have to be made. We suppose that 
a) the dynamical time-scale of the pre-existing 
structure is short compared to the time-scale of 
the infall (i.e. the pre-existing structure is viri- 
alized) and b) that the structure already present 
reacts adiabatically to the infall. In detail, the 
three assumptions about infall of matter at for- 
mation step k are: 

i) Total Energy: The total energy of the in- 
falling mass ruk is conserved. There is no 
energy dissipation and 
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const. 



(6) 



ii) Virialization: The dynamical time-scale of 
the pre-existing structure is short com- 
pared to the time-scale of the infall. The 
infalling mass ruk virializes in the gravita- 
tional potential originating from the pre- 
existing, virialized structure. The corre- 
sponding virial theorem is (e.g. Binney 
&Tremaine 1987) 



2Kk + Wk + Vk 







(7) 



where Kk is the kinetic, Wk the potential 
energy and Vk is related to the external 
gravitational potential. 

iii) Adiabatic Contraction: The pre-existing 
structure reacts adiabatically to the in- 
falling matter deposited within its exten- 
sion. This assumption is based on the fact 



that for periodic orbits / pdq (where p and 
q are then conjugate momentum and coor- 
dinate, respectively) is an adiabatic invari- 
ant (Blumenthal et al. 1986, Flores et al. 
1993). 

To simphfy calculations, the density distribution 
resulting from virialization of the mass rrik is ap- 
proximated by 
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where Rk is the characteristic (virial) radius de- 
termined by (0). 

Before discussing assumption i), i) and iii) in 
more detail, consider the sequence of formation 
steps. The formation of the virialized halo starts 
with the collapse of the mass rape at redshift 
zk- There is no pre-existing structure and the 
collapse according to the three assumptions i), 
ii) and iii) is equivalent to the familiar top-hat 
collapse. The resulting object (structure K) is 
characterized by the density distribution pK @i 
where Rk is determined from (|7|) with Vk = 0. 

The mass mx-i collapses onto structure K 
at redshift zk-i- Its virialized density distribu- 
tion pk-1 is determined by eqs. (0,§) depending 
on pk- The mass enhancement within Rx due 
to the infall of mx-i leads to an adiabatic con- 
traction of the innermost structure K (assump- 
tion iii). It is estimated here by considering the 
quantity RkMk as adiabatic invariant. In prin- 
ciple, this implies a self-similar mass distribution. 
Since we do not attempt to spatially resolve px 
within Rx, this is a reasonable approximation 
here. The contraction of structure K due to the 
infall of rriK-i then is described by 

[mK-i{R'K) + mK]R'K = ^kRk (9) 

where mK-i{R'K) is the fraction of the mass 
rriK-i which is inside the 'contracted' radius R'j^. 
In formation step K — 2, the density distribu- 
tion pk-2 is determined by eqs. (0,^) depending 



now on pK and pK-i- The infall of mass mK-2 
leads to a contraction of structure K — \ and the 
structure K contracts because of mK-2 o,nd the 
contraction of structure K — 1 (which contributes 
to a mass enhancement within Rk as well). 

The formation of the collapsed object is com- 
pleted at z = after the formation steps K — 
3, -fC — 4, . . . , 0. The structure of the virialized 
halo is described as a superposition of concen- 
tric spheres of different but uniform density and 
defined by the set {{nik-, Rk), A; = 0, 1, 2, . . . ,K}. 
The total density profile is 
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Let us stress similarities and differences be- 
tween this model and earlier work. With respect 
to the assumption i) and iii) (no energy dissipa- 
tion, adiabatic response), the model is equivalent 
to the collapse model of Gunn (1977) and models 
(e.g. Hoffman & Shaham 1985, Hoffman 1988) 
based on the stable clustering assumption (Fil- 
more &: Goldreich 1984, Bertschinger 1985). The 
similarity partially extends to the process of viri- 
alization in the sense that a stepwise, separated 
virialization of infalling matter is assumed. How- 
ever there are two main differences: The model 
makes no assumptions about self-similarity and 
it includes the formation history. With respect 
to these differences, the presented model clearly 
adopts a more realistic view of hierarchical struc- 
ture formation. 

But how realistic are assumptions i)-iii) in 
the context of hierarchical structure formation? 
A complete answer to this question lies beyond 
the scope of this letter. The key issue probably 
is assumption i). attempt to semi-analytically 
model heating of the pre-existing structure by in- 
falling matter has recently been made by Nusser 
h Sheth (1998). Their study seems to indicate 
that including energy dissipation (e.g. by drag 
forces, Ostriker &: Turner 1979) does not drasti- 
cally alter density profiles but leads to a slightly 
steeper slope in the innermost region of halos. 



Regarding assumption ii) the time-scale for viri- 
alization which is comparable to the dynami- 
cal time-scale Tmr ~ (47rG/9max)~^'^ is of im- 
portance. As the virialized structure forms and 
the cuspy inner region emerges, r^r rapidly gets 
smaller. From the outcome of the model calcu- 
lation (c.f. Section |3|) one finds that for timestep 
Tin — ^-2 ~ 0.05 and after the first few forma- 
tion steps Tmr ^ Tin- This indicates that it is 
reasonable to assume infalling matter to collapse 
onto an at least pre-virialized structure. Finally, 
it should be stressed that the model uses energy 
considerations to relate mass histories and virial- 
ized density distributions. These considerations 
depend on the gravitational potential which itself 
is much smoother and more symmetric than the 
underlying density distribution (see also Hoffman 
1987). 

3. Results 

The results presented here are based on 50 
Monte Carlo realizations of the formation his- 
tory of halos of mass Mq = lO^^M©, lO^^M© 
and IQi^M© in the context of the SCDM cos- 
mology (O = l,/i = 0.5, (T8 = 0.67). A timestep 
Az = 0.05 has been chosen and the mass limit 
below which the object is considered to be a seed- 
ing object is defined as Mf < 0.05Mo. While the 
time step Az affects the mean mass falling in at 
each formation step, the resulting virialized den- 
sity distributions do not change significantly in 
the tested range Az = 0.01 . . . 0.2. 

3.1. Virialized Density Distribution 

Figure |^ shows a superposition of the virial- 
ized density distributions. Data points mark the 
total density at the various radii R^ ■ The sharply 
dropping density at the outer edge of the distri- 
butions originates from the missing infall of mass 
in excess of the defined halo mass. The radius at 
which this drop occurs roughly corresponds to 
the radius r2oo at which the mean overdensity 
drops below 6 = 200. The radius r2oo is known 



as the characteristic size of a virialized object. 
This implies that for r > r2oo the assumptions 
of the model break down. As a consequence, the 
density distribution at r > r2oo will not be con- 
sidered hereafter. Table || presents the mean of 
r2oo and the corresponding mass M2oo- The scat- 
ter of the innermost data points is due to the 
different formation masses Mf and redshifts Zf 
leading to different initial values of Ri and pi, 
respectively. For radii r < Ri, the structure of 
the density distributions are not resolved. Thus 
the meaningful innermost radius for an individ- 
ual density distribution is r^^^ = Ri. The radial 
range of the model thus is [rrcs,'^20o]- 

Figure g shows the mean virialized density 
profile. Data points mark the mean of the 
50 individual density profiles for the 10^^ Mq, 
10^^ Mq, and 10^^ Mq halos (from left to right). 
Dashed lines indicate the rms scatter. The data 
points mark the radial range [f^es, '^20o]- The best 
fit NFW profile (||) is plotted as solid line. Ta- 
ble [H lists the concentration parameter ro/r2oo- 
Higher halo masses lead to less centrally con- 
centrated halos. The comparison of Figures ^ 
and with Figures 3 and 4 of NFW96 and cor- 
responding concentration parameters illustrates 
the agreement with N-body simulations. 

3.2. Circular Velocities 

Mean circular velocities and the corresponding 
scatter are calculated from Vc{r) = \/GM(r)/r. 
The maximum values for the scatter are found 
roughly at r = r,,,. It is ±13%, ±10% and ±6% 
for the 1O^^M0 lO^^M© and W^^Mq halos, re- 
spectively. These values are in good agreement 
with Eisenstein & Loeb (1995) (see Figure 2 in 
their paper). For the outer regions of the ha- 
los, the model predicts a significantly decreasing 
scatter which reaches values well below 10% at 
1^ = '''200 (c-f. Table |I|). If confirmed by N-body 
simulations, this could be of importance for the 
interpretation of the Tully-Fisher, e.g. in terms 
of its relation to cosmological initial conditions 
(see Eisenstein & Loeb 1995). As a word of cau- 



tion it should be noted that the poor radial den- 
sity resolution at r ^ r^es affects the circular rota- 
tion profiles. Compared to the NFW profile, the 
present model suggests a higher central density 
leading to flatter circular velocity profiles. 

4. Summary and Conclusions 

A simple analytic collapse model has been pro- 
posed for the hierarchical formation scenario of 
virialized cold dark matter halos. As a working 
hypothesis it has been assumed that dark matter 
halos posses a memory of their formation his- 
tory by conserving the total energy of infalling 
matter in each formation step. For the tested 
SCDM cosmology, the virialized mean density 
profiles are in good agreement with numerical re- 
sults. This suggests that energy dissipation does 
not play a major role in the determination of 
the mean shape of virialized density profiles. As 
a consequence, the (nonlinear) virialized density 
distribution is tightly related to the (linear) den- 
sity field of the infall region by the total energy 
of infalling matter. 

An improved model could include more realis- 
tic density distributions pk and probably should 
differentiate between minor and major mergers 
(e.g. Salvador-Sole et al. 1998). While these 
modifications could allow for resolving the radial 
range r < 0.05r200) it is believed that they will 
not drastically alter the results for r > 0.05r2oo- 
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Table 1 
Parameters of the Mean Virialized Density Profiles 
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Fig. 1. — Superposition of the virialized density distributions resulting from the 50 mass histories and 
the collapse model discussed in Section ^. The halo masses are 10^^ Mq, lO^^M© and 10^^ Mq (from left 
to right). 
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Fig. 2. — Mean virialized density distribution (points) and rnis scatter (dashed lines). The sohd hne 
shows the best fitting NFW profile. The outermost data points mark r = r2oo- For a discussion of the 
inner limit c.f. text. 
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Fig. 3. — Scaled mean density profiles (points) and the best fitting NFW profiles of the W^^Mq, W^^Mq 
and W^^Mq mass halos (from top to bottom at r/r2oo < 0.1). For better readability, plotting of the 
rms scatter has been omitted. The scaled profiles agree within one a for r/r2oo > 0.1. Note the different 
behavior of the density distributions for r/r2oo < 0.1. Higher halo masses imply less centrally concentrated 
halos. 



